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1 Introduction 



Understanding the origin of fermion flavor structure, i.e. quark/lepton masses 
and mixing angles, is one of important issues in particle physics. Indeed, var- 
ious types of mechanisms have been proposed to realize realistic mass matri- 
ces of quarks and leptons. Hereafter, we refer to such mechanisms as flavor 
mechanisms. Non-Abelian discrete flavor symmetries are interesting propos- 
als and several types of models with non-Abelian discrete flavor symmetries 
have been constructed pp. 

Origins of such non-Abelian discrete flavor symmetries are not clear in 
4D field theory, but those may be originated from geometrical structures 
of extra dimensional field theories and superstring theories on 6D compact 
spaces. In Ref. [21 [3l H] it has been shown that certain flavor structures 
can be derived from heterotic string models on orbifolds. One of typical 
non-Abelian discrete flavor symmetries, which can appear from heterotic 
orbifold models, is symmetry and it can appear from factorizable orbifolds 
including the Z 2 orbifold. Indeed, several semi-realistic string models with 
the -D4 flavor symmetry have been constructed in Ref. (21 [5], where three 
families consist of D4 singlets and doublets. Thus, it is important to study 
several phenomenological aspects of D4 flavor models. 

Grimus and Lavoura proposed first the -D4 flavor model for the lepton 
sector [B] and subsequently other several models were studied [ZHSJ- Recently, 
the authors proposed a new D4 lepton flavor model [9], which has only a 
single electroweak Higgs field, while the Grimus-Lavoura model as well as 
other models has three electroweak Higgs fields. 

Supersymmetric extension of the standard model (SM) is one of interest- 
ing candidates for new physics beyond a TeV scale. Supersymmetry (SUSY) 
can stabilize the Higgs mass against radiative corrections due to heavy modes 
around high energy scales, e.g. the right-handed Majorana neutrino mass 
scale, the GUT scale, the Planck scale, etc. Supersymmetric standard mod- 
els also have a good candidate for dark matter. Furthermore, in the min- 
imal supersymmetric standard model (MSSM) with a pair of up and down 
Higgs supermultiplets, three gauge couplings are unified at the GUT scale 
Mx — 2 x 10 16 GeV in a good accuracy. 

In supersymmetric models, quarks and leptons have their superpartners, 
i.e. squarks and sleptons. If we have a flavor mechanism to lead to realistic 
masses and mixing angles of quarks and leptons, such a flavor mechanism 
would affect mass matrices of squarks and sleptons. Furthermore, if a spe- 
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cific pattern of sfermion masses is derived by a certain flavor mechanism, 
that would become the prediction of a certain flavor mechanism, which could 
be tested by measuring squark/slepton masses in future experiments. Thus, 
it is important to study patterns of sfermion masses, which are obtained by 
flavor mechanisms leading to realistic fermion masses and mixing angles. Al- 
though squarks and sleptons have not been detected yet, their mass matrices 
are constrained severely by experiments of flavor changing neutral current 
(FCNC) processes [10jo Off-diagonal elements of sfermion mass matrices 
must be suppressed in the super-CKM basis, where fermion masses are diag- 
onalized, when sfermion masses are of order of the weak scale. At any rate, 
it is important to study which patterns of sfermion mass matrices, i.e. soft 
scalar masses and A-terms, are derived from each flavor mechanism. 

In this paper, we study supersymmetric extension of the D4 model of 
Ref. p] as well as the Grimus-Lavoura model. The -D4 model of Ref. [9] 
has a single electroweak Higgs field, while other models have more than one 
electroweak Higgs fields. This difference is important in supersymmetric ex- 
tensions. The former becomes the MSSM with a pair of up and down Higgs 
fields at low energy, but the latter would have more pairs of Higgs fields. 
The latter may violate the gauge coupling unification unless we introduce 
extra colored supermultiplets. That is the reason why we study mainly su- 
persymmetric extension of the D 4 model of Ref. [9] . We evaluate soft SUSY 
breaking terms of sleptons, which are derived in supersymmetric extension 
of the -D4 model [H] as well as the Grimus-Lavoura model. We compute 
soft scalar masses and A-terms for the lepton sector in these supersymmetric 
models within the framework of the gravity mediation of SUSY breaking. 
We examine constraints due to FCNC experiments!! 

The paper is organized as follows. In Section 2, we supersymmetrize 
the -D4 model of [H] and show values of parameters consistent with neutrino 
oscillation experiments. In Section 3, we evaluate soft SUSY breaking terms 
of sleptons, i.e. soft scalar mass matrices and A-terms. We examine FCNC 
constraints on those SUSY breaking terms as mass insertion parameters. 
Section 4 is devoted to conclusion and discussion. In Appendix, we also 
discuss supersymmetric extension of the Grimus-Lavoura model. We evaluate 
soft SUSY breaking terms in the supersymmetric Grimus-Lavoura model and 
show that it is almost the same as the results obtained in Section 3. 



1 See also e.g Ref. [IT] and references therein. 

2 See Ref. [4] for a similar analysis on the quark sector. 
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2 Supersymmetric model with D4 flavor sym- 
metry 

In this section, we study our supersymmetric _D 4 model, based on [9J. The 
D4 symmetry has five irreducible representations, i.e. one doublet 2 and 

four singlets 1 ++ , 1_+ and 1 , where 1 ++ is the trivial singlet and 

the others are non-trivial singlets. Their tensor products are obtained as 

2® 2= 1 ++ © l + _ © 1_ + © lob ® led = le/, (1) 

where a,b,c,d = ± and e = ac and / = bd. 

The three generations of left-handed, right-handed charged lepton and 
right-handed neutrino chiral superfields are denoted by Lj, Rj 1 Nj, (I = 
e, [A, r), respectively, and are assigned to trivial singlets 1++ and doublets 
2. We also introduce a _D 4 doublet, (xi, X2), a D± non-trivial singlet, X-+> 
and a D± trivial singlet, x- They are all SM-gauge singlet chiral superfields. 

Moreover, we introduce additional Z 2 symmetry and we assume that the 
superfields R e , (R^, R T ), x> X-+ have Z 2 -odd charges. The others are Z 2 - 
even. Higgs chiral superfields for both up and down sectors, H u,d , are as- 
sumed to be -D4 trivial singlets and Z 2 -even. These assignments are shown 
in Table 1. Hereafter, we follow the conventional notation that superfields 
and their scalar components are denoted by the same letters. 





L e Lj R e R t N e Nj 


Jju,d 


X X-+ (X11X2) 


z 2 


1++ 2 1 ++ 2 1 ++ 2 

+ + -- + + 


+ 


1++ 1 + 2 

- + 



Table 1: D4 and Z 2 charges. I corresponds to I — /i and r. 
The -D4 x Z 2 invariant superpotential of the charged leptons is given as 

W4 4) = j^L e R e H d X + j^(L»R» + L T R T )H d X + |^( L^R, - L T R T )H d X - + , 

(2) 

up to 5-point couplings, where M p is the Planck scale, M p = 2.4 x 10 18 GeV. 
We assume that scalar components of xi, X2, X an d X-+ develop their vacuum 
expectation values (VEVs) and the -D4 x Z 2 symmetry is broken at a high 
energy scale. Then, the above superpotential (T5]) leads to a diagonal charged 
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lepton mass matrix after the electroweak symmetry breaking by Higgs VEVs, 
v d = (H d ) and v u = (H u ). On the other hand, higher-order operators can 
be also included in the superpotential, 

W L ] = TH L e{R»Xl + RrX2)xH d + • • • , (3) 
p 

and these generate off-diagonal elements of the charged lepton mass matrix. 
As a result, after D 4 x Z 2 symmetry breaking, the effective Yukawa matrix 
{yt)u among the charged leptons and the Higgs field H d is found to be 

(y e a a y ei iOia a - y' efl aa b y ei iOca a + y' efl aa a 

y^aa a - y ! ^ e aa b y^a a - y'^ y^a® 2 + y^a&a 2 ( ■ 
y^ e aa a + y'^ e aa b y^^a 2 - y'^ T a b a 2 y^a a + y'^a b 

(4) 

where a = (xi)/M p , a a = (x)/M p , a b = (x-+)/M p . Here, it has been 
assumed that (xi) = (X2) [9|. The mass matrix of the charged leptons is 
obtained as Mg = v d ye- If ot is allowed sufficiently small compared with 
m e /m T , charged lepton masses are estimated by the diagonal elements of (J4]) 
as follows, 

m e ~ y e a a v d , m M ~ (y^a a - y'^a b )v d , m T ~ (y^a a + y'^a b )v d . (5) 

We need the fine-tuning of parameters, y^, y' , a a and a b , such that we realize 
the mass ratio m M /m r as {y^a a - y'^a b )/{y^a a + y'^a b ) = 0{m^/m T ). That 
is similar to Ref. [B|. Also, the coupling y e must be suppressed to lead to 
ye a a/ '{y 'nOi a + y'u a b) — 0{m e /m r ). Off-diagonal elements of the diagonalizing 
matrix of (jlj) are determined by the parameter a. For example, the (1,2) 
element of the diagonalizing matrix is estimated as 

. (y efl a a -y' a b )av d 
On , (P) 

and it becomes 

6»f 2 ~am r /m M , (7) 

when y efl , y' e ^ ~ 1. At any rate, when a < C(10 -2 ), the mixing angle 6[ 2 is 
smallO Such a small value is not important for the neutrino oscillation in 



3 Small deviations from the diagonal form are important in a certain case (see e.g. 

EE])- 
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our model at the current level of experiments, but important for soft SUSY 
breaking terms as we will discuss. 

Now we consider the neutrino sector. The D4 x ^-invariant superpoten- 
tial is given as 

W® = Vl N e L e H u + V2(N ti L tl + N T L T )H U 

+ y a N e {N, X i + N tX 2) + MiN e N e + M 2 (AyV M + N T N T ), 

up to 4-point couplings. The higher-order operators should be also considered 

W » = W Ne{L » Xl + L ^ RU + W L e( N ^Xl + N tX 2)H U + ■■■. (9) 

With the above superpotential, the Majorana mass matrix, Mr, and the 
Dirac mass matrix, Mq, of the right-handed neutrinos are written down as 

Mi y a M p a y a M p a 
M R = I y a M p a M 2 y b M p a 2 

y a M p a y b M p a 2 M 2 , ^ 

yi yna yna 

M D =V U j 1/21" 1/2 1/23 « 2 

I/21 a y 32 a 2 Vi 

after the -D4 x Z 2 symmetry and the electroweak symmetry are broken. The 
light neutrino mass matrix is obtained as, 



M v = MlM~ l M D =\y z w I , (11) 




where x, y, z and w are given in terms of parameters of Eq. lfroj) . This form 
of My can realize the realistic neutrino mixing without fine-tuning [6j [7J [U [9] . 
The light neutrino masses are written as^ 

'm x \ 
M v = V ( m 2 \V T , (12) 
m-J 



4 These values are obtained at the energy scale Mr and radiative corrections are, in 
general, important to evaluate values at low energy. However, radiative corrections are 



negligible in this D4 flavor model [8]. 



5 



/ cos 9 12 sin6 l 1 2 \ 

V=\-sm9 12 /V2 cos9 12 /V2 1/V2 ), (13) 
• sin 



mi 



m 2 



1 ylr 2 + ylr + y/2yxy 2 y a kr/ (cos 9 i2 sin 12 ) (v u ) 2 



2 r 2 - 2?/ 2 A; 2 r Mi 

1 y 2 r 2 + ylr - \f2y\y2y a krj (cos (9 i2 sin 6> 



12j K) 2 



2 r 2 -2y 2 a k 2 r M 1 ' 

3 r 2 -2y 2 a k 2 r M x ' 1 J 

where we define as r e ^ and fc = ^g^, and neglect higher order terms 
of a such as a 2 M p and a 2 M x appeared in M^ 1 . The justification of this 
approximation and derivation of Eq. (fT4l are given in Ref. The mixing 
angle 6*12 is written as 

cot 2ft, = /f~ Vl , - (15) 

When yi,?/2 and y a are of 0(1), the above mixing angle 9 12 is of 0(1) and 
the effect due to 9\ 2 is negligible in the neutrino oscillation. Thus, the at- 
mospheric neutrino mixing angle is maximal and the Chooz mixing angle is 
vanishing, while the solar neutrino mixing is of 0(1). This also holds if the 
mass matrices were complex, i.e. if the CP violating case would be studied. 

Now, let us consider the realization of experimental values. We use the 
best fit values of mass squared differences and solar mixing angle as [H] 

Am 2 tm = 2.4 x 1(T 3 eV 2 , Am s 2 ol = 7.6 x 1(T 5 eV 2 , 

sin 2 0i2 = 0.32. (16) 

First we consider the simple case with y\ — y 2 — y a — 1- The above experi- 
mental values ffTBl) are obtained by taking^ 

M 1 = 4.9 x 10 15 GeV, M 2 = 6.2 x 10 14 GeV, \a\ = 1.6 x 1CT 3 , (17) 



5 If Mi, M 2 , M p a satisfy the relation Mi +M p a = M 2 , cot26»i 2 = 1/(2^2) is realized at 
this approximation level and the mixing matrix V in Eq. (I13p is the so-called tri-bimaximal 
matrix 1151. 
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for v u ~ 174 GeV0 In this parametrization, the parameter a is sufficiently 
small, so that the charged lepton mass matrix is diagonal as we assumed in 
the beginning of this section. When we vary yi, y 2 and y a around yx, y 2 , y a = 
(9(1), the above experimental f|T6|) values are realized for [9] 

M 2 ~ 0.2 x M u Mj~3x 10 15 x y\ GeV, 

a „ o.OOl x (18) 

Va 

Thus, it is found that the value of a is predicted around O(10 -4 ) — O(10~ 2 ) 
as long as couplings, yi,y 2 ,y a are of 0(1). 



3 Soft SUSY breaking terms in supersym- 
metric D4 model 

We consider soft SUSY breaking terms of slepton mass matrices within the 
framework of supergravity theory, i.e. the gravity mediation^ In our model, 
the -D4 symmetry restricts not only the fermion mass matrices but also the 
scalar matrices. Now we assume chiral superfields to cause SUSY breaking 
by their non-vanishing F-components. The F-components are given by 

F *k = _ e ^ K *k7 ^Qjff + ^W^j , (19) 

where K denotes the Kahler potential, Kjj denotes second derivatives by 
fields, i.e. Kjj = djdjK and K IJ is its inverse. In general, the fields 
in our notation include D4 x Z 2 -singlet moduli fields Z and %i X— h Xi- 
Furthermore VEVs of F<$, k /Q k are estimated as (F$ k /Q k ) = 0(m 3 / 2 ), where 
7723/2 denotes the gravitino mass, which is obtained as 7713/2 = (e K / 2M %W/M%). 

First let us discuss the scalar mass matrices given by using the second- 
order Kahler potential of left-handed and right-handed leptons, 

=a e {<b k )L\L e + cv($ fe )(LtL M + L\L T ) 
+ b e (<S> k )RtR e + b^ k )(RlR fl + RtR T ), 



6 The numerical result has been given in the non-SUSY case [9]. 

7 SUSY breaking may be realized as the gauge mediation or anomaly mediation. These 
mediation mechanisms would lead to the flavor-blind soft SUSY breaking terms. 
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where a 6jAt ($fc), b e ^(§k) are _D 4 x Z 2 -invariant generic functions. Then the 
soft SUSY breaking scalar masses are given by 



m 2 7j K 7j = mi /2 K 7j - \F**\ 2 d 0k d^K 7j + \F**\ 2 d^K lL d* h K m K LM , (21) 

for the scalar fields with the Kahler metric Kjj, where we have assumed the 
vanishing vacuum energy. Evaluating (F® k /<&k) for all fields as (F' s ' k /^j : ) = 
0(1713/2), the slepton mass matrices can be found to be 

/m 2 Ll \ (m 2 R1 \ 

ml = m| 2 , m\= m| 2 , (22) 
V m 2 L2 J V m mj 

where run, m R i = 0(m 3 / 2 ). These contributions give the leading-order of the 
slepton mass matrices. Since the second and third families are -D4 doublets, 
these forms (1221) can be easily expected from the -D4 flavor symmetry and it 
is the prediction of our model that the second and third families of sleptons 
have degenerate masses at this leveljf) 

Eq. (1221) is written in the D4 flavor basis. The super-CKM basis is con- 
venient to examine constraints due to FCNC For example, the (1,2) ele- 
ments are obtained in the super-CKM basis as (^ 2 L ^ R ))u CKM ^ ~ 6[ 2 (^' 2 
m i(i?)2)- T ne mass insertion parameters are defined as 



lm L(/?)i 



(m 2 s(SCKM) 

{°LL(RR))ij - — 2 ' \ Z6 > 

" b SUSY 

where uisusy denotes the average mass of sleptons. The (1,2) elements are 

constrained severely by the /i — >• ej experiments as (b l LL ( RR ))i2 < O(10 -3 ) 
when msusY — 0(100) GeV, while the others have no strong constraints. 
That requires 

e{ 2 < C(10" 3 ). (24) 

For ux = y 2 = y a — 1, we obtain a = O(10~ 3 ) and 6[ 2 = O(10 -2 ), and 
such a parameter region will be ruled out. The value of a is obtained as 
a = 0.001 x ym/Va for y u y 2 , y a ~ 1 and we estimate 0[ 2 = ay ef , for y efM ~ y'^ 



8 When three families consist of a singlet and a doublet under a non-Abelian discrete 
symmetry, a similar structure would appear, that is, two of three families of scalar masses 
corresponding to the doublet would be degenerate at a certain level. See e.g. Ref. [16J. 
However, the exactly same structure as our model has not been predicted in other models. 
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as discussed in Section 2. Possibly, we can get 9{ 2 = (9(10 3 ) in a certain 
parameter region, e.g. y x = y 2 = y ell = y' efl = Vl/a = 0.5. Hence, our 
model is marginal for FCNC constraints, that is, a certain parameter region 
is ruled out already but the remaining parameter region is still wide. Of 
course, if the couplings y efl , y' , which are irrelevant to realization for the 
neutrino oscillation experiments, are suppressed as y e ^y' e ^ ^ ®(X)i a quite 
wider region for yi,y2, y a is allowed. 

In the above estimation we restricted the Kahler potential to K^, but we 
can write more general terms in the Kahler potential after D\ x ^-symmetry 
breaking. The VEVs of Xi are large compared with VEVs of the other fields 
X and X-+ as the result of Section 2. Thus, corrections including Xi are more 
important in correction terms of the Kahler potential. Therefore, we include 
such terms in the Kahler potential. For example, the following terms 



AK (e) = H^±Z L t (L ^t + LrX t) + >^f±Ll {Wxi + LtX2 ) + h . a , (25) 



where ai±(Z) and PLi(Z) are dimensionless generic functions of moduli fields 
Z, appear for the mixing between L e and L^. In addition, we have the 
following corrections 



for the mixing between L M and L T , where ctL±(Z) and PLi(Z) are generic 
dimensionless functions of moduli fields Z. Here, a and cr + _ denote 



and the ellipsis in (126]) denotes terms including other bi-linear combinations 
°f Xii X2 and their conjugates. The contribution of the D4 singlets \i X-+ 
are suppressed by Z 2 symmetry and do not contribute to the leading order of 
off diagonal matrix elements. Also, the Kahler potential for the right handed 
charged leptons R e ,^, T has similar corrections. Including these corrections, 
scalar masses squared are estimated in the D4 flavor basis as 



(26) 



c— = X1X2 - X2X1, 



0+- = X1X2 + X2X1, 



(27) 




m\ x 0(am 2 3/2 ) 
0(aml /2 ) m\ 2 
0(aml /2 ) 0(a 2 ml /2 ) 




(28) 
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m 2 m 0(aml /2 ) 
m\ = | 0(am 2 3/2 ) m 2 R2 0(a 2 m 2 3/2 ) | , (29) 
K 0(aml /2 ) 0{a 2 m\ l2 ) m% 2 

where corrections of 0(a 2 m 3 ^ 2 , alm^ 2 , a%m 3 / 2 , a a abml^ 2 ) are omitted in the 
diagonal elements. Even though we include these corrections, we have almost 
the same constraint due to FCNC as the previous estimation (|24|) . 

Also we examine the mass matrix between the left-handed and the right- 
handed sleptons, which is generated by the so-called A-terms. The A-terms 
are trilinear couplings of two sleptons and one Higgs field, and are obtained 
as [13] 

h IJ D I R J H d = hjjDjRjH 11 - (y £ ) LJ D I RjH d F^K LT d^ k K Tl 

-(y i ) IM D I R J H d F^K MlI d^K Mj , (30) 

where hjj = F® k d$ k (yi)u. Note that effective Yukawa couplings (ye)ij in- 
clude Xi X-+i Xi as Eq. flU). After electroweak symmetry breaking, these 
provide us with the left-right mixing mass squared (rn 2 LR )jj = hjjv d . 

Now, let us discuss the first term in Eq. (1501) . hij. For simplicity, we 
assume that Yukawa couplings are independent of moduli fields Z. For ex- 
ample, the (3,3) element, h 33 is obtained as 

px px-+ 

h ^ = y, Wp + y'^- ( 31 ) 

Then, we can estimate hs3V d = 0{m T mz/2)- Similarly, we obtain the (2,2) 
element as 

px px-+ 

h22 = y »w p - y '^- (32) 

Thus, we evaluate h 2 2V d = (9(m r m 3 / 2 ), because we estimate F x /x, F x ~+ /x-+ 
= 0(7713/2), but F x /x and F x ~+ /x-+ are ; m general, different from each 
other. That may cause a problem. If \hu/yu\ is large compared with slep- 
ton masses, there would be a minimum, where charge is broken [TJjJ^ In 
order to avoid this, we assume that the Kahler metric of x an d X-+ are the 



9 If a decay rate from the realistic minimum to such charge breaking minimum is 
sufficiently small compared with the age of the universe, that might not be a problem. 
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same0 e.g. canonical, and the non-perturbative superpotential leading to 
SUSY breaking does not include x or X-+> i- e - (^x^O = (^x-+^0 = 0- m 
this case, we obtain 

F x/ x = F x- +/x _ + = _ m?j/2i (33) 

and we estimate fi22V d = 0(171^1713/2) ■ Similarly, other elements of hjj are 
estimated as 

/ 0(m e ) 0(m^[ 2 ) 0(m T a)\ 

huv d = m 3/2 0(m^{ 2 ) 0(m M ) 0(m r a 2 ) . (34) 

\C(m r «) C(m r a 2 ) C(m T ) / 

The other terms in Eq. (130]) lead to the same order of A-terms as Eq. fl34|) . 
Thus, we obtain the left-right mixing slepton mass matrix in the D4 flavor 
basis as 

/ 0{m e ) 0(m M 0f 2 ) 0(m T a)\ 
m\ R = h IjV d = m 3/2 0(171X2) 0(m M ) 0(m r a 2 ) . (35) 

\0(m r a) C(m r a 2 ) C(m r ) / 

The pattern of this mass matrix in the super-CKM basis is the same as the 
above. We define the mass insertion parameters for the left-right mixing as 

/_2 s(SCKM) 
fxl \ \" L LR)ij , n\ 

( S Lr)h = -V > ( 36 ) 

''''SUSY 

where (mf jR )( SCKM ^ is the left-right mixing slepton mass squared matrix in 
the super-CKM basis. Only for the (1,2) element of (S l LR )ij, there is a strong 
constraint due to FCNC as (^12 < 0(1(T 6 ) [TO], when m SUSY = 0(100) 
GeV. This constraint also requires 9[ 2 < O(10~ 3 ), which is the same as 
Eq. d2H». 

The soft SUSY breaking terms, which we have studied, are generated at 
a high energy scale such as the Planck scale or the GUT scale. In the above 



10 The D4 flavor structure can be realized in heterotic string models on factorizable 
orbifolds including Z4 [3j. In those heterotic orbifold models, Z?4 non-trivial singlets and 
trivial singlets appear in the same sector and have the same Kahler metric. In such models, 
our assumption would be justified. 
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discussion, we have neglected radiative corrections. The gaugino contribu- 
tions are dominant in radiative corrections to slepton masses, that is, slepton 
masses at the weak scale are obtained by ones at the GUT scale Mx as 



where and are bino and wino masses, respectively. The above esti- 
mation on FCNC constraints does not change drastically when these gaugino 
masses are comparable with slepton masses. If these gaugino masses are quite 
large compared with slepton masses, FCNC constraints would be improved. 

4 Conclusion 

We have studied supersymmetric extension of the D4 x Z2 flavor model of [9] . 
We have evaluated soft SUSY breaking terms about the slepton mass terms. 
It is remarked that the second and third families of slepton masses are almost 
degenerate. The difference is tiny as 0(a 2 rriy 2 , a^rn^ 2 , alm^ 2 , a a abm 2 / 2 ). 
The (1,2) element 9\ 2 of the diagonalizing matrix for the charged lepton 
mass matrix is important for the FCNC constraints, in particular /1 — ► ej 
experiments, although it is not important to realize the neutrino oscillation 
experiments. It is constrained as 6\ 2 < 0(1O -3 ) from the current bound of 
BR(fi — > ej) and our model is marginal. Thus, future improvement on the 
bound of BR(n — > ej), e.g. by the MEG experiment [18] is quite important 
in our model. 

In Appendix, we have also discussed supersymmetric extension of the 
Grimus-Lavoura -D4 x Z 2 flavor model, which leads to almost the same results 
for soft SUSY breaking terms. 

Finally, we give a comment on realization of our SUSY model by string 
model building. The -D4 flavor symmetry can appear from heterotic string 
models on factorizable orbifold models including the Z2 orbifold like Z2 x Zn 
orbifolds [21 El E] , unless one does not introduce Wilson lines, which break 
degeneracy of massless spectra. Indeed, several semi-realistic models have 
been constructed [21 [5] , where three families correspond to D4 trivial singlets 
and doublets. Stringy realization of our flavor structure would be plausible 
from this viewpoint. However, such heterotic orbifold models include only D4 
trivial singlets and doublets, but not Z) 4 non-trivial singlets as fundamental 



m 2 L (Mz) 
m\{M z ) 



m 2 L (Mx) + 0.5M| + 0.04M] 
m 2 R {M x ) + 0.2M], 



(37) 



12 



modes. On the other hand, the D 4 non-trivial singlet x~+ plays an impor- 
tant role in our model. Such a mode could appear as a composite mode. 
Alternatively, -D4 non-trivial singlets as well as trivial singlets and doublets 
can appear as fundamental modes in heterotic string models on factorizable 
orbifolds including the Z4 orbifolds like Z4 x Zn orbifolds. Thus, stringy 
realization on such orbifolds might be alternative possibility. 
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A Appendix 

Here, we discuss the supersymmetric extension of the Grimus-Lavoura model j6]. 
The Grimus-Lavoura non-SUSY model includes three electroweak Higgs fields. 
Thus, we have to introduce three pairs of Higgs superfields for the up and 
down sectors, H^' d (i = 1,2,3). Also we introduce a D 4 doublet {xi,X2)- 
Table 2 shows -D4 and Z2 charges for all fields. 





L e Lj R e Rj N e N T 


TjU,d TjU,d TjU,d 

H 1 ti 2 H 3 




z 2 


1++ 2 1 ++ 2 1 ++ 2 

+ + - + 


1++ 1++ 1+- 

+ + 


2 

+ 



Table 2: D A and Z 2 charges for the supsymmetric Grimus-Lavoura model. / 
corresponds to fi and r. 

Following these assignments, the superpotential which leads the lepton 
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mass matrices is found to be 

ty (3) = y e L e R e H* + y^L^ + D T R T )H* + y'^L T R T - L^H* 

+ {y x L e N e + y 2 (L lt N ft + L T N T ))H? + y 3 (L T N T - L„N^ (38) 
+ y a N e (N^ X i + N tX 2) + M 1 N e N e + M 2 (N^N fl + N T N T ), 

up to 4-point couplings. Furthermore, higher-order superpotential terms like, 

= y e ,( X iR» + X2R T )L e H* + ■ • • , (39) 

should be also considered, because (xi) must be large to lead the realistic 
neutrino mixing. Now the charged lepton mass matrix can be evaluated as 
follows 

/ y e vf M-!/>3)« (Ve^l + yLv$)a\ 
M l = y^vfa Vfi vi - y'vf 0(v d 2 a 2 ) , (40) 

\y„ e vta 0{via 2 ) yA + v'A j 

where (Hf) = vf. If a <C 1 is allowed, eigenvalues of lepton masses can 
be found to be equal to the diagonal elements like Eq. (T4|). We need fine- 
tuning Yukawa couplings and VEVs such that (y^v 2 — y'^vf) / {y^v^ + y'^vf) = 
Oijn^/mr). In addition, we require y e vf to be suppressed compared with 
y^ v 2 + V'a v i to lead to the mass ratio m e /m T , i.e. ye^f/iy^ + v'^ v i) = 
0(m e /m T ). The (1,2) element of diagonalizing matrix 9[ 2 is estimated as 
#12 = {VepVi - y'^vfja/m^ and it reduces to d[ 2 ~ y efl am T /m^ for y efl ~ y' efl 
and v 2 ~ v d . Thus, this is also the same as Eqs. ([6]) and (J7J). 

On the other hand, the Dirac mass matrix, Md, and the Majorana mass 
matrix, Mr, in the neutrino sector are written as 

(Z/12W2 - y'i2^)a (yi 2 v% + y' 12 vf)a 
M D = \ y 21 v\a y 2 v^-y 3 v% 0(v%a 2 ) 

0{v%a 2 ) y 2 v\ + y z vl 

y a M p a y a M p a \ 
M R = I y a M p a M 2 y b M p a 2 , (41) 
y a M p a y b M p a 2 M 2 J 

where (H™) = vf. The above pattern is quite similar to Eq. fflQj) . and in 
particular, the form of Mr is the same as Eq. fflOl) . Thus, similar values 
of parameters lead to realistic results ffl~6l) . i.e., M 1? M 2 = (9(10 15 )GeV and 
a ~ M 2 /M p . Thus the favorable region of a is of O(10~ 4 ) - O(10" 2 ). 




14 



In this model, the soft SUSY breaking terms are also restricted by -D4 x Z 2 
symmetry and expected not to be different from the estimation in Section 3. 
In fact, the scalar mass terms are given by 

jm 2 Ll + 0(a 2 ml /2 ) 0(am 2 3/2 ) 0(am 2 3/2 ) 
m 2 L = 0(aml /2 ) m 2 L2 + 0(a 2 mj /2 ) 0(a 2 m 2 3/2 ) | , (42) 

V 0(aml /2 ) OOa 2 m 2 3/2 ) m 2 L2 + 0{a 2 m 2 3/2 ) 

(m 2 Rl + 0(a 2 m 2 3/2 ) 

m 2 R =\ m 2 R2 + 0{a 2 m\ j2 ) 0{a 2 m\ j2 ) | , (43) 

V 0(a 2 m 2 3/2 ) m 2 R2 + 0(a 2 m 2 3/2 ) 

in the -D4 flavor basis. Thus, we obtain the same constraint on the mass 
insertion parameters due to FCNC as ([21]), i.e. 6[ 2 < 0(1(T 3 ). 

Similarly, the A-terms can also be evaluated. This model would have the 
same problem about the (2,2) element of A-terms as the model in Section 3, 
that is, the (2,2)-element would be of 0(m T m 3 / 2 ) without tuning about the 
coefficient. However, when the Kahler metric of H 2 and H 3 are the same, 
the (2,2) element becomes of 0{m il mj > i2). Then, the left-right mixing slepton 
mass matrix could be estimated as 

/ 0{m e ) 0(6[ 2 m T ) 0{am T ) 
™>lr = m 3/2 0{am e ) 0(m M ) 0{a 2 m T ) | , (44) 
\0{am e ) 0{a 2 m T ) 0{m T ) 

in the -D4 flavor basis, and its form is almost the same in the super-CKM ba- 
sis. Thus, the mass insertion parameter is estimated as (S l LR )i2 = 0(9[ 2 m^/m 3 / 2 ), 
and we have the same constraint as one in Section 3, i.e. 8[ 2 < 0(1O~ 3 ). 

As a result, soft SUSY breaking terms, which are predicted in the super- 
symmetric Grimus-Lavoura model, are almost the same as those obtained in 
Section 3. The difference is the number of Higgs pairs, that is, the super- 
symmetric Grimus-Lavoura model has three pairs of Higgs supermultiplets, 
while the model in Section 2 has only one pair and it becomes the MSSM 
at low energy. The former may violate the gauge coupling unification unless 
one introduces extra colored supermultiplets. 
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